Fluctuations and pattern selection near an Eckhaus instability 
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Qs, \ We study the effect of fluctuations in the vicinity of an Eckhaus instability. The 

; 

classical stability limit, which is defined in the absence of fluctuations, is smeared out 

S 

into a region in which fluctuations and nonlinearities dominate the decay of unstable 
states. The width of this region is shown to grow as D 1 / 2 , where D is the intensity 

O ■ 

of the fluctuations. We find an effective stability boundary that depends on D. A 

•i-H . 

numerical solution of the stochastic Swift-Hohenberg equation in one dimension is 



used to test this prediction and to study pattern selection when the initial unstable 
state lies within the fluctuation dominated region. The asymptotically selected state 
differs from the predictions of previous analyses. Finally, the nonlinear relaxation for 
D > is shown to exhibit a scaling form. 
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Recent experimental studies of Rayleigh-Benard convection in CO2 gas PJ, convection 
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in binary mixtures ||, electro-hydrodynamic convection in nematic liquid crystals ||, and 
Taylor- Couette flow [[|, have detected random fluctuations of thermal origin that are strong 
enough to be analyzed quantitatively. Whereas the effect of thermal fluctuations on ther- 
modynamic equilibrium is now well understood in general terms, much less is known about 
their effect in the vicinity of nonequilibrium instabilities. Pioneering work in this direction 
was done by Graham || , and Swift and Hohenberg || . It was believed, however, that since 
characteristic thermal energies are several orders of magnitude smaller than the character- 
istic driving or dissipative energies involved in many nonequilibrium instabilities, the effect 
of fluctuations would be far too small to be observable. The results of these recent experi- 
ments have renewed interest in the subject since they have opened the possibility of finding 
nontrivial phenomena due to fluctuations close to nonequilibrium instability points. 

We study in this paper the effect of fluctuations on the Eckhaus instability 0. The Eck- 
haus instability is a longitudinal instability often exhibited in systems that display patterns 
which are spatially periodic. It has been exhaustively studied in many systems, including, 
for example, the stability of a set of parallel convective rolls in Rayleigh-Benard convec- 
tion H, the stability of a periodic array of cells during directional solidification [[J, and 
instabilities of stationary standing Faraday waves [|IT]]. We also note particularly detailed 
experimental studies of the Eckhaus instability in electro-hydrodynamic convection in ne- 
matic liquid crystals by Lowe and Gollub ||11|| , and Rasenat et al. ||12|| . Three main issues 



are addressed here: the fluctuation-induced smearing of the classical Eckhaus boundary, 
the nonlinear relaxation of the unstable solution, and the asymptotic periodicity of the new 
stable solution following the instability (this is often referred to as pattern selection). 

In the absence of fluctuations, the Eckhaus boundary of potential or gradient systems is 
located at the line where both the first and second functional derivatives of the appropriate 
Lyapunov functional vanish; hence it separates regions of metastability and instability. This 
is the analog of a spinodal line in a first order phase transition. Our study can be motivated 
by analogy with this latter case: if fluctuations are included in the description of a first order 



phase transition, the spinodal line ceases to exist |13[ . Instead, states that are not stable are 



classified as metastable or unstable with reference to a particular temporal scale of evolution 
(observation time). States that do not decay within such a scale are said to be metastable, 
and unstable otherwise. The often narrow region that separates both types of behavior is 
shifted with respect to the spinodal line, and defines an effective stability boundary known as 
the "cloud point" . We found here related phenomena, namely the smearing of the Eckhaus 
boundary for finite amplitudes of the fluctuations, and the emergence of a different time scale 
for the decay of the unstable state. This scale depends on the intensity of the fluctuations 
and identifies the observable stability boundary. Our analysis of the fluctuation-dominated 
regime is similar in spirit to Binder's determination of the range of validity of a linearized 



description of phase separation Jl4| . He derived the analog of a Ginzburg criterion to describe 
the effect of fluctuations on the long wavelength instability known as spinodal decomposition. 
Finally, we also discuss pattern selection when the initial state is within the transition region. 

Although we expect our analysis to be of wider applicability, we restrict ourselves here 
to the one dimensional stochastic Swift-Hohenberg (sSH) equation |J , 
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which is known to be a model of Rayleigh-Behard convection near onset. The periodicity of 
the stationary solutions of Eq. Ql|) is commensurate with the convective rolls, e is the reduced 
Rayleigh number and £(x,t) is a Gaussian random process satisfying, 

(£(x, t)) = 0, (C(x, t)C{x', f)) = 26(x - x')5(t - f). (2) 

The value of D is proportional to the intensity of the thermal fluctuations in the fluid. For 
D = and e > 0, Eq. ([!]) admits periodic stationary solutions ip q (x) of wavenumber q with 



q-L < q < qi, and q±L = y 1 ± y/e |I5|. However, only those solutions with q_ E < q < q E 
(g_L < q^E < Qe < Ql) are linearly stable. Periodic solutions with q outside this range 
exhibit a modulational instability known as the Eckhaus instability. 

The decay of the unstable solution is triggered by fluctuations. In order to study such 
decay, we linearize Eq. (0) around a stationary periodic solution ip q .(x), with q E < qi < q L . 



Linearization leads to an eigenvalue problem for a spatially periodic operator. According to 
Floquet's theorem, the normalized eigenfunctions f qi ,k( x ) are of the form flEfl , 



oo 

f qi;k ( X ) = y: ^y^+v* , (3) 
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each associated with an eigenvalue X(qi, k), which determines the growth rate of that partic- 
ular eigenfunction. Because the SH equation is real, these eigenvalues satisfy the symmetry 
relation X(qi, —k) = X(qi, k). Also, the invariance of the SH equation under spatial transla- 
tions implies X(qi,k = 0) = 0. Let ±k m be the values of k that maximize X(qi,k) and let 
f m {x) be the real (and normalized) combination of their associated eigenfunctions. Then, 
it is possible to describe the early stages of the decay of the unstable state by a linear 
stochastic equation governing the amplitude u(t) of f m (x), in the eigenfunction expansion 
of if)(x,t) - il> gi (x): 



u{i) = X m u{t) + VDr]{t) . (4) 

r)(t) is a white Gaussian random process which results from the projection of £,{x,t) onto 
f m {x): rj(t) = J dx f m (x)£(x, t), the average of which is zero. When the eigenfunctions in Eq. 
(0) are normalized, the variance of rj turns out to be independent of the initial periodicity 
qi. Equation (|J) is valid until a time t p at which the amplitude u becomes large enough 
such that nonlinearities become important. The solution of Eq. @ with initial condition 
u(t = 0) = is u(t) = h(t)e Xmt , with 

h{t) = VD C ds e- XmS r]{s) . (5) 
Jo 

The amplitude h(t) is a Gaussian random process, averaging to zero, and with a time- 
dependent variance. From its definition it is easy to see that it becomes stationary after a 
time t m ~ A" 1 . After this time, one can replace h(t) by its long time limit, h(t = oo), which 



is a Gaussian variable of standard deviation y (h 2 ) = J D/X m . Within this approximation, 
the calculation of the average time t p at which u(t) crosses a given reference value ut is 



standard and gives [17 



* P = ^ln^^j. (6) 

Two characteristic time scales emerge from this analysis. The time t m signals the end of 
the fluctuation dominated regime and the beginning of a linear deterministic regime. By 
choosing ut equal to a given (small) fraction of the final saturation value of u, t p can be 
interpreted as the time at which nonlinear terms begin to be important. The approximation 
of replacing h(t) by its asymptotic value is obviously true only if t m <C t p , which implies 
u^A m > D. When 

X m « 4 (7) 

there is no clear separation of time scales, and there will not be a distinct linear deterministic 
regime. Nonlinear effects then become important even in the early fluctuation dominated 
regime. 

The arguments give above are of a general nature. Consider now an initially periodic 
state of wavenumber gj > q E . The linear growth rate of the most unstable eigenfunction 
vanishes near the Eckhaus boundary as 

A m ~ (ft - q E ) 2 - (8) 

This can be seen by noting that the symmetry properties of \(qi, k) stated above imply that, 
for small k, X(qi, k) ~ ak 2 + bk A + where a and b are functions of qi and e. b must be 
positive to avoid instabilities with arbitrarily large wavenumber, and a changes sign at the 
Eckhaus boundary qi = q E . Then, to first order, a ~ (qi — q E )- The maximum \ m (<Ii) of 



\(qi, k) is at k m = ±y a/26, which leads to Eq. (f|). In addition, it can be seen that the most 
important wavenumber in the Fourier expansion of the eigenfunction f m (x) is q m = qi — k m . 
From the expression for k m we get 

qi-q m ~ (qi - g^) 3 - (9) 

We now define a transition region of initial wavenumbers qi around q E , within which the 
extent of the linear deterministic regime vanishes. From Eqs. (0) and (|^) we find that this 
transition region is determined by 



< qi - q E < —y/D, (10) 

u T 

where C is a constant of order one. Note that Ut is finite near the Eckhaus boundary, since 
it has to be taken as a finite fraction of the final saturation amplitude of the fastest growing 
mode, and this quantity has no singularity across the Eckhaus boundary. 

These results are restricted to an initially unstable state; hence they only apply to > q E . 
Initial states with < q E which are linearly stable according to a deterministic calculation, 
can also decay because of fluctuations, extending the transition region to qi < q E . The 
smearing of the Eckhaus boundary and the concomitant reduction of the region of linearly 
stable states can be thought of as an effective shift of the Eckhaus boundary to a new value 
q E < q E - It should be stressed, however, that this new effective boundary is only defined 
with reference to a given observation time for the decay of the metastable state. 

We next address the implications of our analysis on the issue of pattern selection. Pattern 
selection refers to the determination of a global wavenumber qj of a configuration obtained 
after a long time, as a function of the initial wavenumber q^. The wavenumber q$ is defined 



here as g/ = vrn/, where n/ is the number of nodes of the configuration per unit length |L8 
If qi lies outside the transition region, the initial stages of the decay of the unstable state 
are well described by linear theory, and the fastest growing mode in the linear regime, q m , 
will mostly determine g/. In this case, the evolution does not lead in observable time scales 
to a periodic configuration of wavenumber q min that minimizes the Lyapunov functional 
associated with Eq. (|T]). In fact, earlier numerical simulations of the sSH equation showed 
that qj w q m when g, is not too close to the Eckhaus boundary |19| . These results were 
consistent with previous simulations of the amplitude equation also in the absence of noise 
20| . On the other hand, if qi lies inside the transition region, nonlinearities and fluctuations 



are likely to alter these conclusions. 

To provide test our predictions we carried out a computer simulation study of the one- 
dimensional sSH equation. Complete details of the algorithm used can be found elsewhere 
| 2T| . We chose e = 0.05625 and discretized Eq. ([!]) on an evenly spaced grid with N = 8192 



6 



nodes and Ax = 27r/32. The equation was integrated forward in time with At = 1 x 10 4 , 
and the results averaged over a number of independent runs, typically 10-20. With this 
choice of Ax, Aq = 2-rr/NAx = 0.004. 

Figure 1 presents the average value of qf as a function of for D/Ax = 0.05, 0.1 and 0.2. 
The dashed line is the value of q m found by numerically calculating the Floquet spectrum 
on the same grid used to solve Eq. (jl|). The analysis also yields qE ~ 1.201. Simulation 
points for qi < qE are also shown in Fig. |l|: although in the stable range according to the 
deterministic calculation, they are seen to decay. Furthermore, for qi far enough from the 
Eckhaus boundary the value of qj obtained from the simulation approaches q m . However 
as qi — qE is decreased, the data deviate from q m . The size of the region where qj deviates 
from q m identifies a transition region that increases with D, in a manner consistent with 
the prediction of Eq. (|10D. By taking the value of Cu^ 1 ~ 0.219, Eq. fllQl) gives a size for 
this transition region extending up to qi ~ 1.25,1.27, and 1.30, for D j Ax = 0.05, 0.1, and 
0.2, respectively, which is consistent with the trend of the data in Fig. [I]. A more detailed 
comparison of Eq. (jT^) with our numerical results has not been attempted since Eq. ( [TDD 
only gives an estimate with undetermined coefficients. 

In contrast with the prediction of linear theory (Eq. (^)), the simulation points for low 
enough q^ lie on straight lines, so that the exponent 1/2 in Eq. (f|) is changed by fluctuations 
to a value close to 1. We have defined an effective stability boundary q~E, a function of D, 
that separates stable from unstable states within our simulation time (t ~ 100, see Fig. 1). 
We define it to be equal to the intersection point between the straight lines in Fig. 1 and the 
line qi = qf. This gives q~E — 1-16, 1.14 and 1.11 for D j Ax = 0.05,0.1 and 0.2 respectively. 
In addition, our data disagree with standard criteria for pattern selection according to which 
1} — Qm or qf = q m i n , when lies in the transition region. 

Lastly we describe our numerical results for the temporal evolution of the dominant 
periodicity of the configuration q(t) after it becomes linearly unstable. The function q(t) 
is defined as q(t) = 7m(t), where n(t) is the number of nodes of ip(x,t) per unit length, 
averaged over independent runs. The decay is found to take place in a characteristic time 

7 



r with an apparent divergence as $ approaches (Je- A similar behavior has been observed 
in simpler zero dimensional stochastic models used to describe the region that separates a 



metastable from an unstable state There, it is possible to characterize the temporal 
evolution in terms of scaling relations. With this motivation in mind, we have found that 
our data is well described by a scaling relation of the form, 

q(t) - QE = (qi ~ q E )f Wi ~ q E ) z ] , (11) 

as shown in Fig. 2 for D/Ax = 0.1. The best scaling is found for z = 1.7. Scaling is also 
observed for other values of D, with z in the range 1.65 - 1.75. 

Detailed comparison of our results with experiments would certainly require at least 
the consideration of a two-dimensional sSH equation, but we expect that the general ideas 
exposed here, namely the existence of a transition region separating stable and unstable 
states, and the identification of an observable stability limit shifted with respect to qE 
and related to an increase of the characteristic relaxation times, are of wider applicability. 
We note, for example, qualitative similarities with some experimental findings. Lowe and 



Gollub [IT] observed the decay of states in the deterministically stable range and their 



evolution towards states of wavenumber intermediate between q m and q m in- in addition, in 



the experiments by Rasenat et al. |L2| , the Eckhaus boundary was determined in a way very 
similar to ours, and it was noted that, in contrast with deterministic theory, the experimental 
results for qj — qi, and the characteristic time scales of the decay of the unstable state as 
a function of q^ were well fitted by straight lines. Further experimental work, especially in 
systems in which the intensity of the thermal fluctuations is relatively large, is certainly 
needed to clarify whether the effects described above can be explained by our theoretical 
arguments. 
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FIGURES 



FIG. 1. 

Final average wavenumber qf as a function of the initial wavenumber for several values 
of the intensity of the fluctuations D/ Ax: O, 0.05; A, 0.1; and, □, 0.2. The straight solid 
lines are linear fits to the data. The dashed line is q m , the wavenumber corresponding to the 
eigenfunction of fastest growth given by linear theory. The error bars in the figure represent 
the 2a confidence interval for the sample studied. 

FIG. 2. 

Scaled average periodicity of the configuration as a function of time after the instability 
(the unsealed values are shown in the inset) for D/ Ax = 0.1. The value of z — 1.7 has been 
used. The values of shown are (from top to bottom in the right of the inset, and the 
corresponding symbols in the figure), q { = 1.119, 1.211, 1.219, 1.230, 1.238 and 1.250. 
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